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We study some of the local CP-odd and magnetic properties of the non-Abelian vacuum with 
use of overlap fermions within the quenched lattice gauge theory. Among these properties are 
the following: inhomogeneous spatial distribution of the topological charge density (chirality 
for massless fermions) in SU(2) gluodynamics (for uncooled gauge configurations the chirality is 
localized on low-dimensional defects with d = 2. .3, while a sequence of cooling steps gives rise to 
four-dimensional instantons and hence a four-dimensional structure of the chirality distribution); 
finite local fluctuations of the chirality growing with the strength of an external magnetic field; 
magnetization and susceptibility of the QCD vacuum in SU(3) theory; magnetic catalysis of the 
chiral symmetry breaking, and the electric conductivity of the QCD vacuum in strong magnetic 
fields. 
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1. Introduction 

Recently both aspects of our studies, properties of the 
QCD vacuum in strong magnetic fields and its non-trivial 
topological structure, attracted much attention in light of on- 
going heavy-ion experiments. The magnitudes of the mag- 
netic fields possibly created there are of order of B ~ m\, 
which may change the critical temperature of the chiral tran- 
sition, several electromagnetic properties of QCD and also 
give rise to new phenomenological effects potentially mea- 
surable experimentally. For the origin and properties of the 
magnetic fields in heavy-ion collisions see the review [p. 

A non-trivial gluonic background may create an imbal- Figure 1: Isosurfaces of the chirality 
ance between numbers of left- and right-handed quarks and distribution, see [||| for details, 
lead to the local strong CP- violation [g, Q]. A typical spatial 

distribution of such imbalance is shown in Fig. |lj where the colors denote some fixed positive (red) 
and negative (blue) values of chirality. It is remarkable, that for uncooled gauge field configurations 
the distribution is irregular and, as we show later, does not follow the instanton pattern, but forms 
low-dimensional structures. 




2. Magnetic-field-induced effects 
2.1 Technical details 

We use a quenched SU (3) lattice gauge theory with the tadpole-improved Luscher-Weisz 
action [gj. To generate statistically independent gauge field configurations we implement the 
Cabibbo-Marinari heat bath algorithm. The lattice size is 14 , and lattice spacing a = 0.105/m. 
All observables we discuss later have a similar structure: (WGW) for VEV of a single quantity 
or ^^2^) for dispersions or correlators. Here G, G\, 62 are some operators in spinor 

and color space. These expectation values can be expressed through the sum over M low-lying 1 
but non-zero eigenvalues iX\ of the chirally invariant Dirac operator D (Neuberger's overlap Dirac 
operator [Q]): 

V where D ¥k = ih¥k, (2-1) 

and 

*<w - e wMy*M-w>m\<m t (2 . 2) 

(iX k + m)(iXp + m) 

where all spinor and color indices are contracted and we drop them for simplicity. The uniform 
magnetic field Fu = S3 = B is introduced as described in [Sj. To perform calculations in the 



'We believe that the IR quantities are insensitive to the UV cutoff realized by selecting some finite number of the 
eieenmodes [Hi 
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Figure 2: Chiral condensate (a) and magnetization/polarization (^Ga^) (b). 



chiral limit we calculated the expression ( |2.1[ ) or ( |2.2| ) for a non-zero m and averaged it over all 
configurations of the gauge fields. Then we repeated the procedure for other quark masses m and 
extrapolated the VEV to m — > limit. Simulation details for the Sections || and |] are similar, but 
for two colors, and are described in the mentioned there papers. 



2.2 Chiral condensate 

In this section we present our results for the chiral condensate 

E = -(0|W|0), 



1), 



(2.3) 



as a function of the magnetic field B. The result is shown in Fig. 2(a) . The general tendency for £ 
to grow with B has been obtained in various models (see []9|] for a review) and usually referred to 
the magnetic catalysis of the chiral symmetry breaking 2 . We perform the fit of the results by the 
following function 



L fit (B) 



, eB 



(2.4) 



where En = £(0)- The obtained fitting parameters are 

£ = [(228 ±3) MeV] 3 , A B = (1.31 ±0.04) GeV, v = 1.57±0.23. (2.5) 

The value of the condensate in absence of the magnetic field £(0) and the exponent v have reason- 
able values, compare with e.g. [12, 13]. 

2.3 Magnetization and susceptibility 

In this section we calculate the quantity 



(2.6) 



2 A non-monotonic behavior in the vicinity of T c has been recently observed [ jlo[ flX| ] and named as inverse magnetic 
catalysis. 
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Figure 3: Fluctuations of the chirality (a) and electromagnetic current/charge (b). 

where o a p = — \y a , yp] and %{F) is a coefficient of proportionality (susceptibility), which depends 
on the field strength. 



This quantity was introduced in [14] and can be used to estimate the spin polarization of 



the quarks in external magnetic field. The magnetization can be described by the dimensionless 
quantity \l= %■ qB, so that 



(2.7) 



The expectation value ( |2.6| ) can be calculated on the lattice by ( |2.1| ) with & = <J a g . The result 



is shown in Fig. |2(b)| . We can see, that the 12-component grows linearly with the field, which 



agrees with []14J. This allows us to find the chiral susceptibility %(0) = Xo from the slope of the 
curve. After a linear approximation (Tcr^^) = &"eB, where 3 



we obtain Of* = (172.3 ±0.5) MeV and 



fit 



-4.24±0.18GeV~ 2 . 



(2.8) 



(2.9) 



This value fits well into the range of present theoretical estimations (see [15] for a review and 



[16] for recent results). We also have to mention a well known analytic result obtained by the 
OPE combined with the idea of pion dominance [17], supported by two holographic ones [ fl8| , |T|], 
but giving a too high value comparing to our results. This disagreement seems to be an important 
puzzle to solve, because it will lead to a better understanding of the pion dominance assumption 
and the large N c limit. 



2.4 Evidences of the chiral magnetic effect 

One example of a new effect mentioned in the Introduction is the chiral magnetic effect (CME), 
which generates an electric current along the magnetic field in the presence of a nontrivial gluonic 



in our simulation we calculate the magnetization of the d-quark condensate, thus q - 
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Figure 4: Conductivity at various temperatures (a) and at small temperatures and various lattices (b). Here 
(°xc> <%) is the conductivity along (transverse to) the magnetic field. 

background [Q, |J. This effect may naturally take place in heavy-ion collisions and is at the moment 



under active experimental search [20, 21, 22] (see also a review [El] on the interpretation of the 



experimental data). Lattice evidences of the effect can be found in [24, 26, £5], 27]. Here we 
implement the procedure from [|4j for the SU (3) case and study the local chirality 

(2.10) 



ps(x) = y(x)y 5 V(x) = p L {x) - p R {x) 
and the electromagnetic current 

j ll (x)=^(x) Y ^(x). 



(2.11) 



The expectation value of the first quantity can be computed by ( |2.l[ ) with G = 75 and with G = 7^ 
for the second quantity. Both VEV's are zero, as expected, but the corresponding fluctuations 
obtained from (22) are finite and grow with the field strength (see Fig. [|]). Here we use the "IR" 
subscript to emphasize, that we subtract from the quantity its value at B = 0: 

(Y) IR (B) = 1^ d 4 x(Y(x)) B -^J d 4 x(Y(x)) B=0 (2.12) 

V V 

One can inteipret the enhancement of the current fluctuations either as short-living quantum fluc- 
tuations or as a charge flow. We measured the conductivity of the vacuum to argue in support of 
the latter. 



3. Conductivity of the QCD vacuum 

Electric conductivity can be extracted from the correlator of two vector currents, 

Gij{x) = J J 3 f(j,-(0,0)7 ; -(^T)> 



(3.1) 



Following [28], let us define the spectral function p (w) which corresponds to the correlator (3.1) 



dw r,r \ 1 \ -i w \ w cosh(w(r- ~U) 
-K( W ,,) PU ( W ), wtth K M = - S U(*) 



(3.2) 
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Figure 5: Chiral IPR (defined in the text) as a function of lattice spacing (a), and the dimensionality of the 
chirality distribution at various cooling stages (b). 



where T is the temperature. The Kubo formula for the electric conductivity then reads [ ]29| , 

Pij (g>) 



lim 

£0^0 



AT 



(3.3) 



In the limit of the weak time-independent electric field Ek, one has (jj) = OiuEk- We performed 
the SU(2) quenched lattice computations (for details see [^]) and obtained results shown in Fig. [| 
One can conclude from the results, that the QCD vacuum is an isotropic conductor in the de- 
confinement phase, and the conductivity does not depend on the strength of the magnetic field 
(a = 15 ± 2MeV for T = 350 MeV > T c ). In the confinement phase the conductivity in the trans- 
verse directions is zero within error range, while it is a growing function of B -field in the longi- 
tudinal direction. In other words, the strong magnetic field at small temperatures turns the QCD 
vacuum from an insulator to an anisotropic conductor. 

Such a behavior may also lead to the enhancement in the soft photon and dilepton production 



rates in the direction transverse to the magnetic field [ 25 ] 



4. Fractal dimension of the chirality distribution 



We studied some properties of the spatial distribution of chirality ( |2.1C ), such that the volume 
occupied by the distribution and the fractal dimension of that volume, for various lattice spacings 
and cooling stages. To measure the former we used the inverse participation ratio (IPR), 



IPRg = N 



5>( 



Via 



V distribution 



(4.1) 



A=0 



where N is the total number of sites of the lattice and the brackets [...]a=o denote an averaging over 
all zero modes and further averaging over all gauge field configurations. In general, the IPR is equal 
to the inverse fraction of sites occupied by the support of a distribution. The fractal (Hausdorff) 
dimension d can be extracted from the IPRs at various lattice spacings a, 



c 



(4.2) 
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Lattice parameters and simulation details can be found in The results are shown in Fig. |5| Our 
conclusion is that the chirality distribution (and hence fermionic zero modes) has zero volume in 
the continuum limit (since IPR grows with decreasing a) and its fractal dimension equals d = 2.3 



for uncooled configurations (in agreement with [|30j |31j [32], [330), tending to d = 4 after the cooling. 



This behavior supports the idea of [J34[], stating that the Yang-Mills vacuum structure is different 
for different resolutions of the measuring procedure. The low dimensional structure of the vacuum, 
if true beyond the probe quark limit, may lead to a new phenomenology relevant for the heavy-ion 
experiments [35,|36|, [37], 38]. 
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